Curvas Epidemicas con incidence2 y ggplot2

Jorge Mario Estrada Alvarez

Usaremos en esta clase la base del paquete outbreaks

Por favor descarga e instala el el paquete y extraigamos la base denominada: fluH7N9_china_2013

# nombrala de la misma manera para que podamos seguir la clase
h7n9 <- outbreaks::fluH7N9_china_2013

Dos formas de producir Epicurvas

  1. Mediante el paquete incidence2 (descarga e instala)

  2. Mediante ggplot2 de manera mas personalizada

Siempre antes de TODO asegurese de:

  1. La variable fecha que se utilizara para la epicurva debe tener formato fecha

  2. Asegurese que tiene un rango de fechas acorde al brote enun rango esperado

Mediante el paquete incidence2

Se requieren 2 pasos para trazar una curva epidémica con el paquete incidence2:

  1. Crear un objeto de incidencia (utilizando la función incidence())
  • Proporcionar los datos
  • Especificar la columna de fecha a date_index =
  • Especificar el ´interval =` en el que deben agregarse los casos (diario, semanal, mensual..)
  • Especificar cualquier columna de agrupación (por ejemplo, sexo, hospital, resultado)
  1. Representar el objeto de incidencia
  • Especificar etiquetas, colores, títulos, etc.
library(incidence2)

con un dataset individual (linelist)

epicurva <- incidence(x,
  date_index,
  groups = NULL,
  counts = NULL,
  count_names_to = "count_variable",
  count_values_to = "count",
  date_names_to = "date_index",
  rm_na_dates = TRUE,
  interval = NULL)

https://www.reconverse.org/incidence2/articles/incidence2.html

plot(x = ,      
  width = 1,
  colour_palette = vibrant,
  border_colour = NA,
  na_colour = "grey",
  alpha = 0.7,
  fill = NULL,
  legend = c("right", "left", "bottom", "top", "none"),
  title = NULL,
  angle = 0,
  size = NULL,
  nrow = NULL,
  n_breaks = 6L,
  show_cases = FALSE)

Para datos agrupados:

  date_of_onset outcome n
1    2013-02-19   Death 1
2    2013-02-27   Death 1
3    2013-03-07   Death 1
4    2013-03-08    <NA> 1
5    2013-03-09   Death 1
6    2013-03-13 Recover 1
incidence(x = ,
          date_index = "date_of_onset",
          counts = "n",
          groups = "outcome")